(N 



cn 

o 



US-FT/7-01 
June, 2001 



The contribution of off-shell gluons 

to the structure functions F2 and F£ 
and the unintegrated gluon distributions 



A.V. Kotikov 

o_ 

O ■ Bogoliuhov Laboratory of Theoretical Physics 

Joint Institute for Nuclear Research 
^ \ 14-1980 Dubna, Russia 



(N ; A.V. Lipatov 

Department of Physics 
^ ■ M.V. Lomonosov Moscow State University 

^n '. 119899 Moscow, Russia 



G. Parente 



Departamento de Fisica de Particulas 
Universidade de Santiago de Compostela 
1^ ' 15706 Santiago de Compostela, Spain 

p,' N.P. Zotov 

(D ■ 
^ ' D. V. Skobeltsyn Institute of Nuclear Physics 

M. V. Lomonosov Moscow State University 

<^ '. 119899 Moscow, Russia 

Abstract 

We calculate the perturbative parts of the structure functions F2 and F^ for a 
gluon target having nonzero transverse momentum squared at order a^. The 
results of the double convolution (with respect to the Bjorken variable xb 
and the transverse momentum ) of the perturbative part and the unintegrated 
gluon densities are compared with HERA experimental data for Fg^. The 
contribution from F^ structure function ranges 10 -=- 30% of that of F2 at the 
kinematical range of HERA experiments. 



PACS number(s): 13.60.Hb, 12.38.Bx, 13.15.Dk 



1 Introduction 

Recently there have been important new data on the charm structure function (SF) 
F|, of the proton from the HI [|I|, |] and ZEUS g |[ Collaborations at HERA, which 
have probed the small-xe region down to xb = 8 x 10"'' and xb = 2 x 10"'', respectively. 
At these values of Xb, the charm contribution to the total proton SF, Ff, is found to 
be around 25%, which is a considerably larger fraction than that found by the European 
Muon Collaboration at CERN |^ at larger xb, where it was only ~ 1% of Ff'. Extensive 
theoretical analyses in recent years have generally served to confirm that the Fg data can 
be described through perturbative generation of charm within QCD (see, for example, 
the review in Ref. [§ and references therein). 

In the framework of DGLAP dynamics [^ §] there are two basic methods to study 
heavy flavour physics. One of them [^ is based on the massless evolution of parton distri- 
butions and the other ||10| on the boson-gluon fusion process. There are also interpolating 



schemes (see Ref. [jTT| and references therein). The present HERA data [^ ^, ^ ^ for the 
charm SF F| are in good agreement with the predictions from Ref. pi| . 

We note, however, that perhaps more relevant analyses of the HERA data, where the 
Xb values are quite small, are those based on BFKL dynamics [|12| (see discussions in the 
review of Ref. |T3| and references therein), because the leading ln{l/xB) contributions 



are summed. The basic dynamical quantity in BFKL approach is the unintegrated gluon 
distribution ipg^XjkjJ {fg is the (integrated) gluon distribution multiplied by xb and k± 
is the transverse momentum) 



cQ' dk^ 



/" dk 
—^ (Pg{xB, k\) (hereafter g^ = -Q^, k'^ = -k\ 



± ;> 



which satisfies the BFKL equation. 
We define the Bjorken variables 

Xb = QV(2w) and x = Q^/{2kq), (2) 

for lepton-hadron and lepton-parton scattering, respectively, where p^ and A;^ are the 
hadron and the gluon 4-momentums, respectively, and q^ is the photon 4-momentum. 

Notice that the integral is divergent at the lower limit and so it leads to the necessity 
to consider the difference fg{xB,Q'^) — fg{xB,Qo) with some nonzero Ql (see discussions 



in Sect. 3), i.e. 

fg{xB,Q^) = fg{xB,Ql)+ -rr'fai^B^kl) (3) 

In our analysis below we will not use the Sudakov decomposition, which is sometimes 
quite convenient in high-energy calculations. However, it is useful to have relations be- 
tween our calculations and the results, where the Sudakov decomposition has been used. 
The corresponding analysis will be done in the next Section. Here we only note that the 
property fc^ = —k^ (see Eq. (P) comes from the fact that the Bjorken parton variable x 
in the standard and in the Sudakov approaches coincide. 

Then, in the BFKL approach the SFs F|^(xb,<5^) ^ire driven at small xb by gluons 
and are related in the following way to the unintegrated distribution (pg{x, fc^.)- 

'■1 dx r dk'j 

Ixg X _ 

The functions C|^(x, Q^, m'^, kj_) may be regarded as the structure functions of the off- 
shell gluons with virtuality k'j_ (hereafter we call them as coefficient functions) . They are 
described by the quark box (and crossed box) diagram contribution to the photon-gluon 
interaction (see Fig. 1). 

The purpose of the article is to calculate these coefficient functions Clj^{x, Q^, rnl, k\) 
and to analyze experimental data for F2{xb,Q'^) by applying Eq. (^ with different sets 
of unintegrated gluon densities (see Ref. |l^) and to give predictions for the longitudinal 
SF Fl{xB,Q'). 

It is instructive to note that the diagrams shown in Fig. 1. are similar to those of 
the photon-photon scattering process. The corresponding QED contributions have been 



FldxE, Q') = t - [^ Cidx, Q\ ml kl) Vg{xB/x, kl), (4) 

Jxa X J K \ 



calculated many years ago in Ref. [|T5| (see also the beautiful review in Ref. |T6[)- Our 
results have been calculated independently and they are in full agreement with Ref. [^ 
(see Appendix B). However, we hope that our formulas which are given in a more simple 
form could be useful for others. 

The structure of this article is as follows: in Sect. 2 we present the basic formalism 



of our approach with a brief review of the calculational steps (based on Ref. |T^). The 



connection of our analysis with the Sudakov-like approach is also given. Later, we present 





Figure 1: The diagrams contributing to T^^, for a gluon target. They should be multiphed 
by a factor of 2 because of the opposite direction of the fermion loop. The diagram (a) 
should be also doubled because of crossing symmetry. 



the results for the two most important polarization matrices for off-shell gluons into the 
proton. In Sect. 3 and 4 we give the predictions for the structure functions F2 and F£ 



for two cases of unintegrated gluon distribution functions (see Ref. [|T4|) used, which are 
shortly reviewed. In Appendix A we show the basic technique for the evaluation of the 
required Feynman diagrams. Appendix B contains the review of QED results from Refs. 
flEj , . In Appendix C we consider the limiting cases, when the values of the quark mass 
or the gluon momentum are equal to zero and also when the value of the photon "mass" 
Q^ goes to zero. 



2 Approach 

The hadron part of the deep inelastic (DIS) spin-average lepton-hadron cross section can 
be represented in the form [|: 

Ff,u = e^^{q) Fl{xb,Q'^) +d^^{q,p) F2{xb,Q'^), (5) 

where g^ and p^ are the photon and hadron momenta, 



e^i^l?) = fl'M^ - -^ and d^^{q,p) 



g^^, + 2xb—^ 5 — + Pt^P,.^^ 

qz g^ 



and Fk^XB, Q^) (hereafter k = 2,L) are structure functions. 

The tensor F^^, is connected via the optical theorem with the amplitude of elastic for- 
ward scattering of a photon on a hadron T^i,{q,p), which may be decomposed in invariant 
amplitudes Tf:{xB,Q^) by analogy with Eq. (|^). 

Let us expand the invariant amplitudes in inverse powers of xb'- 

Tk = E - ^M (6) 

n=0 ^-^B/ 

The coefficients T^ „ coincide (for even n) with the moments M^ „ of the SF F^: 

n,n = Mk,n= fdzz'^-^Fkiz^Q^) (7) 

Jo 

2.1 Evaluation of coefficient functions 

We would like to note that the previous formalism can be replicated at parton level 
by replacing the hadron momentum p^^ by the gluon one k'^ and the Bjorken variable 
Xb by the corresponding x. Then, the hadron part of the deep inelastic spin-average 
lepton-parton cross section can be represented in the form 

F',u = e,M) m^, Q') + d,uiq, k) Ff (x, Q^), (8) 

where F^{x, Q^) are the structure functions of lepton-parton DIS. 

As in our analysis we only consider gluons, the unintegrated gluon distribution into 
the parton (i.e. into the gluon) v^^(x, A;^.) should have the form 

^i{x,ki) = 6{i~xmki), 



^Hereafter we consider only one-photon exchange approximation. 



where ^{k\) is a function of k^. 

The parton SF -F|'(x, Q"^) and the amphtudes T|'(x, Q"^) at the parton level obey equa- 
tions similar to Eq. (|^) and Eq. (|^) with the replacement xb — ^ x. Then, they are 
connected via optical theorem: 

nix^Q") = T.[-) dz z--'Faz,Q') {n = 2m) (9) 

^o\xJ Jo 

Thus, the coefficient functions Cf (x, Q'^,k'^) of the parton SF F^i^x, Q"^) 

r dk"^ 
FU^.Q') = J ^ Cl^{y,Q\mlkl)e{xo-x) 0{kl), (10) 

can be obtained directly using the amplitudes T^(a;, Q^) at the parton level 

r dk"^ 
TUx^Q') = J ^Cl,{y,Q^mlkl)0{kl), (11) 

in the following way: 

C',{x,Q\e) = Y.[-) dzz-~'Ci{z,Q\e)Q{z,-z), (12) 

where we have extracted a kinematical factor Q{zq — z). 

As it was already discussed we will work with gluon part alone, keeping nonzero val- 
ues of quark masses and the gluon virtuality. The corresponding Feynman diagrams are 
displayed on Fig. 1. The coefficient functions Cf (x, Q^, A;^) do not depend on the target 
type. So, it can be calculated in photon-parton DIS and used later in the photon-hadron 
reaction (see Eq. (|)). 

2.2 Connection with the Sudakov-Uke approach 

One of the basic ingredients in the Sudakov-like approach is the introduction of an 
additional light-cone momentum n^ with n^ = and {np) = 1. 
The gluon momentum k^^ can be represented as 

k"^ 4- k'^ 
k^ = iv''+ I, ^ n^ + fc^ (13) 

with the following properties 

p^ = 71^ = [pkx) = {nkr) = 0, (np) = 1 (14) 



where the four-vector k!^ contains only the transverse part of k'^, kx = (0, A;_l,0), i.e. 
k"^ = —k\. and ^ = Xb/x is the fraction of the proton momentum carried by the gluon 
(see Eq. (§. 

To study the relations between the "usual" approach used here and the Sudakov-like 
one, it is convenient to introduce the following parametrization for the vector n^ (see Ref. 



0) 



2t 
n^ = -^{xEP^ + qn (15) 

It is easy to check that the properties in Eq. ( [I^ ) are fulfilled. 

Then, for the scalar product {kq) we have in the Sudakov-like approach: 

k^ + k"^ 
(kq) = ^ (pq) + — — ^ (nq) 



Q 



2 



2x 
li k'^ = —k']_, then it follows from Eq. ([T6| ) 



1- X 



Q' 



(16) 



X 



2{kqy 
that agrees with Eq. (Q). Also from Eq. (|T^) it follows 

k^' = xp" + fc^, (17) 

where x is the fraction of the proton momentum carried by the gluon. 

2.3 Feynman-gauge gluon polarization 

As a first approximation we consider gluons having polarization tensor (hereafter the 
indices a and (3 are connected with gluons and yU and v are connected with photons)^: 

pa/3 _ _gap ^^g^) 

This polarization tensor corresponds to the case when gluons do not interact. In some 

sense the case of polarization is equal to the standard DIS suggestions about parton 

^In principle, we can use here more general cases of polarization tensor (for example, that one based 
on the Landau or unitary gauge). The difference between them and Eq. (n^ is ~ k°' and/or ^ k^ and, 
hence, it leads to zero contributions because the Feynman diagrams in Fig.l are gauge invariant. 



properties, excepting their off-shell property. The polarization in Eq. ( ]TB| ) gives the main 
contribution to the polarization tensor we are interested in (see below) 



pa/3 
^BFKL 



rv \ rv \ 



X 



,py^ 



(19) 



ki -k^' 

which conies from the high energy (or kx) factorization prescription []19|, ^, ^ 0. 

Contracting the photon projectors (connected with photon indices of diagrams in 
Fig.l.) 

P!.'J = -l9,. and P(^)=4x^^ 

with the hadronic tensor F^iy, we obtain the following relations at the parton level (i.e. 
for off-shell gluons having momentum k^) 



P' Clix) 
/3^ Clix) 



)C 



K 



fW I jL_ f (2) 

2/32 
4bx'f^'^ + il±3^ /(2) 



(20) 
/C /(2) + 46a;2/^2 Cf , (21) 



where the normalization factor K, = el as{Q'^)/{4:n) x, 

/C/«,z = l,2 






and /J^ 



Ahx' 



-k^ IQ'^ = k\IQ'^ > 0, a = rr? IQ'^ . The kinematical factor z^ 



which appear in Eq. (12) is 



Applying the projectors P'^} to the Feynman diagrams displayed in Fig.l, we obtain 
the following results for the contributions to expressions 



-2/? 



1 - 2x(l + 6 - 2a) [1 - x(l + 6 + 2a)] /i 



+ (2a-6)(l-2a)xV2 



(23) 



/ 



(2) 



\xi5 



(1 - (1 + h)x) - 2x\hx{\ - (1 + h)x){\ + 6 - 2a) + afiA /i 



+ 6a;^(l -(l + 6)x)(2a-6) /2 



(24) 



■^We would like to note that the BFKL polarization tensor is a particular case of so-called nonsense 
polarization of the particles in t-channel makes the main contributions for cross sections in s-channel at 
s ^ oo (see, for example, Ref. |22] and references therein). The limit s ^ oo corresponds to the small 
values of Bjorken variable xs, that is just the range of our study. 

^ The contributions of individual scalar components of the diagrams of Fig.l (which come after eval- 
uation of traces of 7-matrices) are given in Appendix A. 



where 

Aax 



andg 



^ (l-(l + 6)x) 



1 , 1 + /3/3 



/i = ^t; ^^z ;7^' f- 



pp 1-pp l-P^P 



2R2 



The important regimes: k"^ = 0, m'^ = and Q^ ^ are considered in Appendix C 
The Q^ = hmit is given in Sect. 2.5. 

2.4 BFKL-like gluon polarization 



Now we take into account the BFKL gluon polarization given in Eq. (|T9|). As we 
already noted in previous subsection, in these calculations we did not use Sudakov de- 
composition and, hence, the hadron momentum p" is not so convenient variable in our 
case. Thus, we represent the projector Pbfkl ^^ ^ combination of projectors constructed 
by the momenta k°' and g". 

We can represent the tensor F"^ in the general form: 

F"^ = Ag""^ + Bq''q'^ + Ck''k'^ + D[k''q'^ + q''k^), (25) 

where A, B, C and D are some scalar functions of the variables y, a and b. 

iFiom the gauge invariance of the vector current: k^'F"'^ = k^F'^^ = we have the 
following relations 

Ce = -[A + D [kq]] , B {kq) = - D k^ (26) 

If we apply the BFKL-like projector Pbfkl ^-nd use the light-cone properties given in 
Eq. (p^), we get the simple relation 

PfFKLF"-^ = D (kq) (27) 

The standard projectors g°''^ and q°'q'^ lead to the relations 

g'^ppc.p ^ sA + D^(3^ (28) 



(29) 



/?2 

q<^qPpaP = f- 



SA + D^P^ 



2x 



^ We use the variables as defined in Ref. 



/^From Eqs. 



^), (^) and (^ we have 



:ia/3 



^P' 



,a/3 



and the BFKL-hke projector Pbfkl can be represented as 



BFKL 



1 1 

2;^4 L' 



^2^a/3 _ ^26x 



2gV 



Q^ 



(30) 



(31) 



In the previous section we have already calculated the contributions to coefficient 
functions using the first term within the brackets in the r.h.s. of Eq. (^). Repeating 
the above calculations with the projector ~ g"g^, we obtain the total contribution to the 
coefficient functions which can be represented as the following shift in the results given 
in Eqs. 






^2,BFKL\^) 



where 



/ 



(1) 



/ 



-/? 



(2) 



f(l) 
J BFKL 



f(2) 
J BFKL 



1 

J' 






2/(2) 



^L,BFKL\^)'^ 



P'f 






(32) 



1 - a;(l + b) 



X 



2 x(l-x(l + 6))(l + 6-2a) + a/3M /i 



x(l-x(l + 6))(l-2a)/2 



/ 



(2) 



4/3 (l-(l + 6)x)^ 



2 - (1 + 26x2) /i - bx^ h 



(33) 



(34) 



For the important regimes when /c^ = 0, 



ra 



and Q2 — > 0, the analyses are given 



in Appendix C. 

Notice that our results in Eqs. (p3D and (0) should coincide with the integral repre- 
sentations of Refs. |19|, |25[ (at Q^ -^ there is full agreement (see following subsection) 



with the formulae of Refs. [|T9|, ^ for photoproduction of heavy quarks). Our results in 
Eqs. (0) and (|3^) should also agree with those in Ref. |2^ but the direct comparison is 



quite difficult because the authors of Ref. |Q used a different (and quite complicated) 
way to obtain their results and the structure of their results is quite cumbersome (see 
Appendix A in Ref. ||2^). We have found numerical agreement in the case of F2{xb, Q"^) 
(see Sect. 3 and Fig. 4). 



2.5 Q^ = limit and Catani-Ciafaloni-Hautmann approach 

We introduce the new variables s, p and A which are useful in the limit Q^ —* 0: 

s = — , p = Aax = — ^, A = ox = „„ X 



X 



Q' 



(35) 



and express our formulae above as functions of p and A at small x asymptotic (i.e. small 

When X = we have got the following relations: 
for the intermediate functions 



P' 



1, P' 



Y^ = /3^ h = hn^-^^m, /2 = ^(l-A)(36) 



/ 



(1) 



2P 



1 + p - ^\ L0) - {1- p) + (2 + p -2L 



A 



+ 2{l0)-1)a' 



xf 



(1) 



-2/3 



2(1- A) -pL(/3) 



/ 



(2) 



xf 



(2) 



(37) 
(38) 



and, thus, for the coefficient functions 



CI 



0, 



c. 



'^ J BFKL-i 



(39) 



where 



Jbfkl 



2/9 



1 + P-^)l(/3)-(1 + p) + (8 + p-(2 + 3p)L(/3))a 



+ 2(l(/3)-4)a2 



(40) 



We note that the results coincide exactly with those from Catani-Ciafaloni-Hautmann 
work in Ref. [|l^, |2^ (see Eq. (2.2) in Ref. |2^) in the case of photoproduction of heavy 
quarks. The 0{x) contribution in the Q^ -^ limit is given in Appendix C (see subsection 
C.3). 
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3 Comparison with F2 experimental data 

With the help of the results obtained in the previous Section we have analyzed HERA 
data for SF F| from ZEUS § and HI [§ collaborations. 



3.1 Unintegrated gluon distribution 

In this paper we consider two different parametrizations for the unintegrated gluon 
distribution [|14|. Firstly, we use the parametrization based on the numerical solution of 
the BFKL evolution equation [^ (RS-parametrization). The solution has the following 



form 26 



^(x,k 



Oi 



02 



C„ 



as + «4 
as 



^0 I I Qo 



as + X 



«2 + «3^+|^ 
1/2 , 



+ a X + 



P 



k'^ ' \k^ I e + ln(l/x) 

1 - aex^nn (itVag) (1 + ana;)(l - x)"«+"^o'"('^'/"«\ (41) 



where 



a 



(42) 



f 1, if k"^ < go(3^)> 

1 qo{x)/k\ if k^ > ql{x) 

The parameters (ai — an, a, (3 and e) were found (see Ref. |^|) by minimization of the 

differences between the l.h.s and the r.h.s. of the BFKL-type equation for the unintegrated 

gluon distribution $(x, k'^) with 'Qq =4 GeV^. 

Secondly, we also use the results of a BFKL-like parametrization of the unintegrated 

gluon distribution $(x, A;^.)/^^)) according to the prescription given in Ref. [^. The 



proposed method lies upon a straightforward perturbative solution of the BFKL equation 
where the collinear gluon density xG{x,fi'^) from the standard GRV set |28| is used as 



the boundary condition in the integral form of Eq. (1). Technically, the unintegrated 
gluon density is calculated as a convolution of the collinear gluon density G(x, /x^) with 
universal weight factors 



$(x, kl, /x^) = /' g{v, kl ^?) - G(-, ^?) d^ 

Jx n n 



where 



g{r],kl,fi'^ 



a^ f Jo(2y'a,ln(l/r/)ln(^VA;i)), if kl < fi 
Vk± \/o(2J«,ln(l/r7)ln(A;i//x2)), if kl > fi 



(43) 



(44) 
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20xun01 where Jq and Iq stand for Bessel functions (of real and imaginary arguments, re- 
spectively), and as = Stts/vr. The parameter «« is connected with the Pomeron trajectory 
intercept: Ap = as41n2 in the LO and Ap = as4:ln2 — Na"^ in the NLO approximations. 



where iV is a number, A^ ~ 18 [^-[0. However, some resummation procedures proposed 
in the last years lead to positive value of Ap ~ 0.2 — 0.3 (see Refs. |3^, ^ and references 
therein). 

Therefore, in our calculations with Eq. (^3|) we only used the solution of the LO 
BFKL equation and considered Ap as a free parameter varying it from 0.166 to 0.53. This 
approach was used for the description of the pr spectrum of D* meson electroproduction 
at HERA [p3 where the value for the Pomeron intercept parameter Ap = 0.35 was 



obtained 0. We used this value of Ap in our present calculations with /i^ = Qg = 1 — 4 
GeV^. 

3.2 Numerical results 

For the calculation of the SF F2 we use Eq. (H) in the following form: 

F^{x,Q') = [' ^Cl{-,Q',0)yG{y,Ql) + 
Jx{i+4a) y y 

+ E /f ■ V /if" *i ciC-, Q^ *i) «(!/, *i, c&l (45) 

Here $(j/, k\, Ql) = -p-{pg{y, k\, Ql) and the changes xb —^ x,x —^ y were done in com- 
parison with Eq. (^. The coefficient function Cf (^) is given in Eq. ( pO]) with fsp^L ^^^ 
fsFKL instead of f^^^ and f^'^\ respectively (see Eq. (^)). The functions f^^^ and f^"^^ 
in Eq. ( pOD are given by Eqs. (|^) and (p^ , respectively, and the functions f^^^ and /*^^^ 
can be found in Eqs. (^) and (0), respectively. 

The integration limits in Eq. (^) have the following values: 



.(1) _ ^n , A^ , 91 

^2 I ' ymax 



yi;,L = x(l + 4a + ^), y«, = 2x(l + 2a); 



^'^'^ - Ql kf2^. = (f-(l + 4a))Q2; 



Lmm 



(2) _ 9^/1 , o^N „(2) 



max 



Vmin = 2x(l + 2a), y, 

1,2(2) _ ^2 ,2(2) _ ^2. /.^N 



^ Close values for the parameter Ap were obtained, rather, in very different papers (see, for example, 
Ref. [pSJ) and in the L3 experiment p6{ . 
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• ZEUS 1996 -1997 



0.6 



O 



0.^ 



0.2 



o.e 



Q^=1,8GeV^ 



i 

1 



* 1,2,3 



Q^=4 GeV^ 




Q^=7 GeV^ 



1,2,3 




0.^ 



0.2 



0.6 



0.^ 



0.2 



0.0 



Q^=11 GeV^ 



Q^=18GeV^ 





Q^=30 GeV^ 




^ ^ Q^=60GeV^ 




^ ^^ Q^=130 GeV^ 




10'^ 10"^ 10'^ 10'^ 10° 



10"^ 10"" 10"^ 10"^ 10"'' 10° lO'" 10'^ 10'^ 10'^ 10° 



Figure 2: The structure function F2{x, Q"^) as a function of x for different values of Q^ 
compared to ZEUS data |3]. Curves 1, 2 and 3 correspond to SF obtained in the standard 



parton model with the GRV [^ gluon density at the leading order approximation and 
to SF obtained in the kr factorization approach with RS ||2^ and BFKL (at Qq = 4 
Ge"V^) [^ parametrizations of unintegrated gluon distribution. 
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a A 



HI Preliminary 



or 



0,3 



0,2 



0,1 



0,8 



0,6 



a A 



0,2 



Q^= 1,5 GeV 




Q^= 3,5 GeV 




Q^= 6,5 GeV 



1,2, 3 




0,0 



Q^= 12 GeV 



Q^= 25 GeV 



1 2 





Q^=60 GeV 




-5^-3-2-1 0^-3-2-1 ^ -3 -2 -1 

10 ar Ig -i' Ig '^^ 



Figure 3: The structure function ^2^(0;, Q^) as a function of x for different values of Q^ 
compared to HI data 0]. Curves 1, 2 and 3 are as in Fig. 2. 



The ranges of integration correspond to the requirement of positive values in the 

I) and (|3^) and also obey to the 



arguments of the square roots in Eqs. (|5), O), (0y 
kinematical restriction {z = x/y) < zq with zq from Eq. 

In Figs. 2 and 3 we show the SF F| as a function x for different values of Q^ in 
comparison with ZEUS [Q and HI |^ experimental data. For comparison we present the 
results of the calculation with two different parametrizations for the unintegrated gluon 
distribution $(x, k\, Ql) in the forms given by Eq. (||) and Eq. (||) at Qo = 4 GeV^. 

The differences observed between the curves 2 and 3 are due to the different behaviour 
of the unintegrated gluon distribution as function x and k^. 

We see that at large Q^ {Q"^ > 10 GeV^) the SF F2 obtained in the fcr factorization 
approach is higher than the SF obtained in the standard parton model with the GRV 
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Figure 4: The structure function F2{x, Q"^) as a function of x at Q^ = 60 GeV^ compared 
to ZEUS data ||. Curves 2 and 3 correspond to RS (at Ql = A GeV^) [||] and BFKL 
(at Qo = 1 GeV^) ||2^ parametrizations obtained with our off mass shell matrix and ones 



from Ref. 24 



gluon density at the LO approximation (see curve 1) and has a more rapid growth in 
comparison with the standard parton model results, especially at Q^ ~ 130 GeV^ |37|. 
At Q^ < 10 GeV^ the predictions from perturbative QCD (in GRV approach) and those 
based on the kx factorization approach are very similar [] and show the disagreement with 
data below Q^ = 7 GeV^ 0. Unfortunately the available experimental data do not permit 
yet to distinguish the kx factorization effects from those due to boundary conditions ||2^ . 
Fig. 4 shows the structure function F2 at Q^ = 60 GeV^ obtained with two different 
gluon densities, i.e. RS (at Qq = 4 GeV^) and BFKL (at Qq = 1 GeV^) parametrizations. 
The difference between curves 2 and 3 are mainly due to the different Ql value used (as 
we have already shown in Figs. 2 and 3, the difference due to the parametrizations is 

^This fact is due to the quite large value oi Qq = A GeV^ chosen here. 

*A similar disagreement with data at Q^ < 2 GeV^ has been observed for the complete structure 
function F2 (see, for example, the discussion in Ref. |3^ and reference therein). We note that the 
insertion of higher-twist corrections in the framework of usual perturbative QCD improves the agreement 
with data (see Ref. Bfl]) at quite low values of Q^. 
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Figure 5: The structure function F2{x, Q'^) as a function of Q^ for different values of x 
compared to ZEUS data M. Curves 1, 2 and 3 are as in Fig. 2. 
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Figure 6: The structure function -F2^(x, Q"^) as a function of Q^ for different values of x 
compared to HI data 0]. Curves 1, 2 and 3 are as in Fig. 2. 

essentially smaller). From Figs. 2, 3 and 4, we note that the difference between the 
kT factorization results and those from perturbative QCD increases when we change the 
value of Qo in Eq- (3) from 4 GeV^ to 1 GeV^ ||3^. In addition, for each case presented in 
Fig. 4. we have done the calculations with our off mass shell matrix elements and those 
from Ref. [^] |. The predictions are very similar and cannot be distinguished on curves 
2 and 3. 

For completeness, in Figs. 5 and 6 we present the SF F2 as a function Q^ for different 
values of x in comparison with ZEUS and HI |^ experimental data. 

^We would like to note that Ref. |gj| contains several slips: the propagators in Eq. (A.l) and the 
products (pq) in Eqs. (A. 4) and (A. 5) should be in the denominator, the indices 2 and L in Eqs. (A. 4) 
and (A. 5) should be transposed. 
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Predictions for F^ 



To calculate the SF Ff^ we have used Eq. p5| ) with the replacement of the coefficient 
function Cf by C£, which is given by Eq. (^I]) . 

In Fig. 7 we show the predictions for Ff^ obtained with different unintegrated gluon 
distributions. The difference between the results obtained in perturbative QCD and from 
the krp factorization approach is quite similar to the F^ case discussed above. 

The ratio W^ = FI/F2 is shown in Fig. 8. We see R^ ~ 0.1 -^ 0.3 in a wide region 
of Q^. The estimation of R'^ is very close to the results for R = Fl/{F2 — Fl) ratio (see 
Refs. 01-@1)- We would like to note that these values of i?^ contradict the estimation 



obtained in Refs. 0, Q. The effect of R'^ on the corresponding differential cross-section 
should be considered in the extraction of F2 from future more precise measurements. 

For the ratio R'^ we found quite flat x-behavior at low x in the low Q^ region (see Fig. 
8), where approaches based on perturbative QCD and on kr factorization give similar pre- 
dictions (see Fig 2, 3, 5, 6 and 7). It is in agreement with the corresponding behaviour of 
the ratio R = Fl/{F2-Fl) (see Ref. [|0|) at quite large values of Ap (Ap > 0.2-0.3). 
The low X rise of R'^ at high Q^ disagrees with early calculations [Q in the framework 
of perturbative QCD. It could be due to the small x resummation, which is important 
at high Q^ (see Fig 2, 3, 5, 6 and 7). We plan to study in future this effect on R in the 
framework of kx factorization. 



5 Conclusions 

We have performed the calculation of the perturbative parts for the structure functions 
F2 and F£ for a gluon target having nonzero momentum squared, in the process of photon- 
gluon fusion. The results have quite compact form for both: the Feynman gauge and a 
nonsense (or BFKL-like) gluon polarizations. 

We have applied the results in the framework of kr factorization approach to the 

analysis of present data for the charm contribution to F2 (-Fg ) and we have given the 
^"At small values of Ap, i.e. when x^^^ ~ Const, the ratio R tends to zero at x ^ (see Ref. p4|). 
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Figure 7: The structure function Fl{x,Q'^) as a function of x for different values of Q^. 
Curves 1, 2 and 3 are as in Fig. 2. 
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predictions for F£. The analysis has been performed with several parametrizations of 
unintegrated gluon distributions (RS and BFKL) for comparison. We have found good 
agreement of our results, obtained with RS and BFKL parametrizations of unintegrated 
gluons distributions aX Qq = A GeV^, with experimental F2 HERA data, except at low 
Q^ {Q^ < 7 GeV^) 13. We have also obtained quite large contribution of the SF F£ at low 
X and high Q2 (q2 > gg GeV2). 

We would like to note the good agreement between our results for F2 and the ones 
obtained in Ref. [|^ by Monte-Carlo studies. Moreover, we have also good agreement 



with fits of HI and ZEUS data for F2 (see recent reviews in Ref. [Q and references 



therein) based on perturbative QCD calculations at NLO. But unlike to these fits, our 
analysis uses universal unintegrated gluon distribution, which gives in the simplest way 
the main contribution to the cross-section in the high-energy limit. 

It could be also very useful to evaluate the complete F2 itself and the derivatives of F2 
with respect to the logarithms of 1/x and Q^ with our expressions using the unintegrated 
gluons. We are considering to present this work and also the predictions for F^ in a 
forthcoming article. 

The consideration of the SF F2 in the framework of the leading-twist approximation 
of perturbative QCD (i.e. for "pure" perturbative QCD) leads to very good agreement 
(see Ref. [^ and references therein) with HERA data at low x and Q^ > 1.5 GeV^. 
The agreement improves at lower Q^ when higher twist terms are taken into account p^ 



As it has been studied in Refs. ||3^, ^, the SF F2 at low Q'^ is sensitive to the small-x 
behavior of quark distributions. Thus, the analysis of F2 in a broader Q^ range should 
require the incorporation of parametrizations for unintegrated quark densities, introduced 
recently (see Ref. [^ and references therein). 

The study of the complete SF F^ should also be very interesting. The structure 
function Fl depends strongly on the gluon distribution (see, for example, Ref. @), 
which in turn is determined [^ by the derivative dF2/dlnQ'^ . Thus, in the framework of 
perturbative QCD at low x the relation between F^, F2 and dF2/dlnQ'^ could be violated 

by non-perturbative contributions, which are expected to be important in the Fl case (see 

^^It must be noted that the cross section of inelastic cc— and 66— pair photoproduction at HERA are 
described by the BFKL parametrization at a smaller value of Qq (Qq = 1 GeV^) ||l4| . 
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Ref. |]SD|). The application of present analysis to Fl will give a "non pure" perturbative 



QCD predictions for the structure function that should be compared with data |41, 4^ 
and with the "pure" perturbative results of Ref. PU 
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A Appendix 

Here we present the contribution to the amplitude of the DIS process from scalar diagrams 
Q in the elastic forward scattering of a photon on a parton. In analogy to Eq. (||) one 
can represent any one-loop diagram of the elastic forward scattering 

Dn^,m,m,m, = [ 7^^ d"'^ (h) dT^k - k,) d^' (q + k,) d'^^q + k, ~ k), (A.l) 

•' [2-K) 

where d{k) = {k"^ — m^)~^, in the form 1^: 

■^ /1\"~ /•l/(l+4a+b) ^__^ 

n=0 ^3;/ ^0 

For application of Eqs. ( |A.1| ) and ( |A.2|) to F2 and Fl coefficient functions, only even n 
are needed. We choose k{n) so that k{2n) = 1. 



Below we rewrite Eq. (|A.2|) in the symbolic form: 



Then, we can represent the needed formulae by: 



-'^^These diagrams appear after calculation of the traces of diagrams in Fig.l. 

^•^This method is very similar to that in Refs. |^, |5^ in the case of zero quark masses. Usually the 
consideration of nonzero masses into Feynman integrals complicates strongly the analysis and requires 
the use of special techniques (see, for example, Ref. [Q) to evaluate the diagrams. Here it is not the 
case, the nonzero quark masses only modifies the upper limit of the integral with respect to the Bjorken 
variable (see the r.h.s. of Eq. ( |A.2| )). 
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. ine loops: 

-Doiio — 


-^ k{n) = 1, 


foiw{z) = 1 


-Diooi ~ 


-^ k{n) = {-ir, 


/iooi(2) = 1 



(A.4) 



2. The triangles: 

-Diiio 
-Dioii 



^0111 — ^ ki^n) = 1. /iiio(^) = /oiii(2;) = -z fi{z) 

Diwi — Hn) = {-ir, fiouiz) = hwiiz) = -z fiiz) 



(A.5) 



3. The boxes: 



-D2110 
Dun 



k{n) 
k{n) 



;i + (-!)' 



hiw{z) = -Az^ f2{z) 
fnniz)=Az'Mz) 



(A.6) 

(A.7) 



We would like to note that the result for the second box -Dim is very similar to ones 
for the triangles. In the case of massless quarks this property has been observed in Refs. 



n, n 



B Appendix 

We compare the results obtained in Sect. 2 and Appendix A with well known formulae 



obtained in earlier works (see Refs. [|T5| , p!6|). 
Following Ref. 



let us consider the kinematics of virtual 77 forward scattering. 
According to the optical theorem (see Sect. 1) the quantity F^^""^ is the absorptive part 
of the 77 forward amplitude, connected with the cross section in the usual way. (The 
expression of the amplitude in terms of the electromagnetic currents is given in Sect. 1). 
In the expansion of F^'^'^^ into invariant functions one should take into account Lorentz 
invariance, T-invariance (symmetry in the substitution fiu <-h> ajS) and gauge invariance 
as well, i.e. Q 



Qi 



Tp^uaP __ V Tp^uaP __ a rp^iuafS __ rp^iuafS 



iB.l] 



The tensors in which F^^'^^ is expanded can be constructed in terms of vectors gf , 
^2 and the tensor g^^ . In order to take into account explicitly gauge invariance, it is 
convenient to use their linear combinations: 



Q? = 


-t' 


Q2 Qi 2 ' 
L qf J 


R>"' = 


= R''^' = 


- g^" + X~^ 



Q". 







qr 



{ (11(12] 

ql 



^2uv 



{qiq2){qiq2 + qM - qiq2q2 - ?29i?i 



(B.2) 



(B.3) 



^Sometimes we replace <? — > gi and fc — ^ (72 with the purpose of keeping the symmetry 1 ^^ 2 in our 
formulae in the first part of Appendix B. 
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where 



^ = (^1^2)^ -qI-qI 



(B.4) 



The unit vectors Qf are orthogonal to the vectors gf and the symmetrical tensor R^^ 
is orthogonal both to gf and ^2; i-6- to Qi and Q2: 



q>^Q^l = ql^Qi^ = 0, q^R'^'' = Q'^R'^" = 0, 



(B.5) 



We note, that R'^" is a metric tensor of a subspace which is orthogonal to gf and 
g2- In the c.m.s. of the photons, only two components of R'^'^ are different from 

(i?^^ = Ryy = 1). 

The choice of independent tensors in which the expansion is carried out, has a high 
degree of arbitrariness. We make this choice so that these tensors are orthogonal to each 
other, and the invariant functions have a simple physical interpretation: 



-iM na n^ /n/3 



+ QtQ''^QlQ'2Fss + 



R^'^R'"^ + R^'^R"'' - R^'^'R"^ 



TT 



R^''QtQ2 + R^^QiQ2 + {^^^ ^ a/?) 






(B.6) 



+ 



^/.^^a/3 _ ^^,p^v 






R'^'Q^QI - R^^QIQ'^ + {^^v ^ aP) 



Tpa 



The dimensionless invariant functions Fab defined here only depend on the invariants 
^^ = (<?! + Q'2)^, <?! and ^2- The first four functions are expressed through the cross 
sections aat {a,b = S,T for scalar and transverse photons, respectively). The amplitudes 
F^f^ correspond to transitions with spin-flip for each of the photons with total helicity 
conservation. The last two amplitudes are antisymmetric. 



We would like to represent the results of Ref. [|15| in terms of our functions, introduced 
in Section 2. 

First of all, we return to the variables introduced in Sect. 2. Then, we have 



X 



4x2^' ^' 



46x2 



^^ + Wx'^' 



Q'2 






^' + 2^^' 



The results of Ref. |]T5| have the form 



15. 



I3F, 



TT 



Ax(3 



1 + 4(a - 1 - 6)T + 126T^ - U - 8aV + 2(2a -1-b 



PF, 



TS 



+2bx{l -x{l + b)))T + 2Ab'T' /i + 6x^2 



-16x/? 



T - 2x<^ ax - 6(1 + 2x(a - 1 - b))T - Q¥T^ \fi 



[B.7) 



"'^^The original results of Ref. |15 contain an additional factor [— Tra^/Q^] in comparison with Eq. 
(B.8), that has to do with the different normalization used in our article (see Eq. ([l|)) and in Ref. p6[. 
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+bx^{6a-b + 6bT)Tf2 



pFsT = -I6xp 



T - 2x1 ax - (1 + x{2a - 1 - b))T - 6b^T^\fi 



+x^i6a-b + 6bT)Tf2 



PFss 

PFts 
PFtt 

ms 



UbxT^l3 

8xp 



2-{l + 2bx^)h-bx^h 

2 



2aT + (1 - a^)^ + GfeT^ + x^|2(l + a + b)-b- 2b{2a - I - b)T 

P 

-12b'T'^f\^, 



2x-3T + x'^\2a-l-b + 66t|/i 



4/3 



pF^T = - lefe^/^x^T/? 



x-4:T+\2T-x\fi-bx^f: 



2/i - /2 



where 



Doing the needed projections on Eqs. ( |B.6| ) we can express the above functions as 
combinations of /(i), f^^\ /(^^ and /(2) (see Sect. 2). 



x(l -x{l + b)) 



^'Fss 



/(^), /3^F^^ = /3V« + ^/(^) 



/3^F, 



T5 



/(2) _ jm 



P'F, 



TT 



^2/(1) _ ^(1) 



(B.8) 



The coefficient functions calculated in Sect. 2 can be expressed as combinations of 
Fab iA,B = S,T). 

For non interacting gluons: 



P^C2 = /C 



P'^Cl = /c 



Fss + FsT — 2{Ftt + Fts) 



Fss - 2Fts + 46x2 (f^t - 2Ftt 



For the BFKL projector: 



/3 C2,BFKL — A^ 



/5 Cl,BFKL — ^ 



Fss + -^5T + -Prs + Ftt 



Fss + ^T5 + 46x2 (f^t + Ftt 



(B.9) 



(B.IO) 
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C Appendix 



Here we consider the particular cases: fc^ = 0, m'^ = and Q^ 
compare with others. 



which are relevant to 



C.l The case /c^ = 

When fc2 = 

C|(x) = /C 
with 



f'' + ^ /(^) 



and Ci{x) = /C/(^) 



(C.l) 



/ 



(1) 



-2(3 



HP) 



(l - 2a;(l - a;)(l - 2a)] - (l - 2x(l - 2a) + 2x^(1 - 4a^ 



/(2) = 8x(3 



where 



[l-x)- 2xa L{(3) 



P'' = l 



Aax 



(l-x) 

and the function L{(3) has been defined in Eq. (^). 

Equations (|C.1|) -( |C7^ ) coincide with the results of Ref. |^ 
Indeed, we have 



(C.2) 



(C.3) 



C| = /C(-2)/3 



l-4x(2-a)(l-x) 



1 - 2x(l - 2a) + 2x^(1 - 6a - Aa^) L{j3) 



Ci = ICSxp 



1 - x) - 2xa L(/3) 



(C.4) 
(C.5) 



The consideration of the BFKL projector does not change the results given above 
because the additional terms (see Eq. (|3T|)) are proportional to k"^ and they are negligible. 
The expression in Eq. ( |U.5| ) also coincides with the corresponding result in Ref. |]25| (see 
Eqs. (A17,A18)). 

C.2 The case m^ = 

When m? = the coefficient functions Cf (x) are defined through /*^^^ and Z*^^-* (see Eqs. 
(PPP, (pip) being in this case 



/ 



(1) 



2 - 1 - 2x(l + b)+ 2x^(1 + 6)2 L(/3) 



(C.6) 
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/ 
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For the coefficient functions themselves, we have 
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(C.^ 



(C.9) 



In the case of the BFKL projector, the coefficient functions Cl{x) are defined by Eqs. 
(PI), (0) and dg) with the replacement f^^ -^ ffp^i as in Eq. (|2|). In Eq. (|2|) the 
expressions for /*^*^ can be found in Eqs. (|C.6|) , (|C.7|) while for J*^*^ are given by: 



and, thus. 
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For the coefficient functions C^ ^pj^j^i^x) we have the following results: 
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(C.13) 



(C.14) 
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C.3 The case Q^ ^ 

Using the definitions in Eq. (p5l), when x -^ we have got the following relations (at 
0(x)): 

for the intermediate functions: 
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(C.19) 
(C.20) 



f(i) 



Similarly to Eq. (|C.18|) , the coefficient functions in Eq. (|39D and the functions fsFKL 
in Eq. ( ^0]) have the additional terms proportional x. 
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1 - A)(2 - 3A) - (1 - A)2 - zA 10) 



(C.23) 



(C.24) 
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